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Magnetospheric Birefringence Induces Polarization 
Signatures in Neutron-Star Spectra 
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ABSTRACT 

We study the propagation of polarization light through the magnetosphere of neu- 
tron stars. At intermediate frequencies (the optical through the infrared), both the 
birefringence induced by the plasma and by quantumelectrodynamics influence the ob- 
servation polarization of radiation from the surface of the neutron star. Because these 
two processes compete in this regime, we find that polarization observations can con- 
strain the properties of the neutron-star magnetosphere, specifically the total charge 
density. We calculate both the phase-resolved and the phase-averaged polarization 
signatures induced by magnetospheric birefringence. 
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1 INTRODUCTION 

The propagation of polarized radio waves away from neu- 
tron stars (NS s) has been well understood for many years 
(Barnard 1986). At this energy scale, it is reasonable to ne- 
glect the coupling of the light to vacuum-based QED ef- 
fects. Lai and Ho (2003ab) have discussed the mechanisms 
by which polarized X-rays are affected by the atmospheres of 
NSs. X-ray photons couple the greatest in the atmosphere as 
the plasma densities and field strengths are sufficiently high. 
Lai and Ho included vacuum-based effects in their analysis. 
Similar effects exist at IR/optical wavelengths when propa- 
gation of the waves through the magnetosphere (possessing 
a lower-density plasma and lower strength magnetic field) 
is considered. This paper presents a model and results that 
show that there is a wavelength-dependent polarization sig- 
nature in the IR/optical wavelength as a result of interac- 
tions between radiation and the magnetosphere. 

For most radio pulsars emission from the magneto- 
sphere itself rather than the surface is thought to domi- 
nate the optical and near infrared emission. However, sev- 
eral radio pulsars exhibit optical radiation from the surface 
{e.g. PSR B0656+14 , PSR B0950+08 and Geminga) (e.g. 
IZharikov et al.l 2004). Radio pulsars are not the only mem- 
bers of the neutron-star menagerie. Tbe optical and near- 
infrared emission from the radio quiet neutron stars such 
as RX J1856. 5-3754 and RX J0720.4-3125(e.#. iPons et alJ 
120021: iMotch et al.l l20031 most likely comes from the sur- 
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face, and possibly the emission from anomalous x-ray pul- 
sars emer ges from a surface layer heated b y magnetospheric 
currents (Th ompson fc Beloborodovl 120041) . With these ob- 
jects in mind let us proceed. 

To make progress we will make several simplifying as- 
sumptions. First, we assume that the charge carriers both 
positive and negative have the mass of an electron and that 
the plasma is cold. A pair plasma is necessarily hot. Second, 
we assume that the Goldreich- Julian model is correct. The 
Goldreich- Julian picture describes the equilibrium charge 
density for a rotating neutron star with its magnetic pole 
aligned (or anti-aligned) with the spin axis. Here we will 
assume that the two axes are not necessarily aligned (es- 
pecially in § 13.21 . Finally, the Goldreich- Julian equi librium 
even in its realm of applicability is unstable (see iMichell 
2004). In spite of these deficiencies, the Goldreich- Julian 
model provides a well-understood starting point to demon- 
strate that observations of polarization of the surface radi- 
ation from neutron stars may provide a window to diagnose 
the properties of the plasma and to encourage further re- 
search to probe more realistic models. 



CALCULATIONS 
EVOLUTION 



POLARIZATION 



iKubo fc Naeatal il983f) outline how the polarization of light 
is affected by a birefringent and dichroic medium: 



-=!lx s + (T x s) x s. 

dz 



(1) 



This equation is valid for a medium with dispersion, emis- 
sion and absorption but not scattering. The polarization is 
described in terms of normalized Stokes parameters s = 
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(si,s 2 ,s 3 ) = (Q/I, U/I, V/I). The vectors fi and T are the 
birefringent and dichroic vectors. We have extended the ex- 
pressions of iHevl fc Shavivl (2000) for a permeable medium 
to include the effects of the plasma: 

((Ve£b — ceVb) sin 2 9cos2<j> \ 
-(Ve£b - £B»7s)sin 2 6 sin 20 , (2) 
S cos0(e s + i(77B-e B )sm 2 6>) J 

where u> is the angular frequency of the radiation; 9 is the 
angle between the direction of propagation and the mag- 
netic field; and cf> is the angle between the projection of the 
magnetic field onto the plane perpendicular to the direction 
of propagation, and the component of the magnetic field in 
the direction of the magnetic moment of the NS. t/b and t B 
are components of the permeability tensor, and t/e and e E 
are components of the dielectric tensor, which are, to first 
order in the fine-structure constant (a/): 

l + ±*-(4X, (3) 
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for B <C Bq and hu> <C m e c 2 where v = ui^/u 2 , and it 
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eB/(m e c) is the cyclotron frequency and 
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(9) 



is the plasma frequency under the assumptions mentioned 
in the introduction. 

To get the second equality we have as- 
sume the Goldreich- Julian (1969) plasma density, 
ncj = (QB cos 9')/(2nec). Q = 2n/P is the spin fre- 
quency of the star and 9' is the angle between the rotation 
axis and the local magnetic field. B is the strength of the 
magnetic field perpendicular to propagation direction and 
Bq w 4.4 x 10 13 G is the critical QED field. In the range 
of NSs modeled, it was found that the dichroic vector 
had little effect on the polarization, and was neglected for 
subsequent results presented here. 

A fourth-order Runge-Kutta integrator with variable 
step size was used to determine the evolution of the po- 
larization vector from the surface of the star to the light 
cylinder, along a Schwarzschild null geodesic. The polar- 
ization vector s was initially set to be in the direction of 
the plasma-dominated birefringent vector, as the initial po- 
larization modes are determined by the plasma-rich atmo- 
sphere of the NS. 

We set the initial conditions by specifying the inclina- 
tion angle a (not to be confused with the fine structure con- 
stant otf) between the optical axis and the magnetic axis, 
and the angle /3 between the plane containing the optical 
and magnetic axes and the plane containing the photon tra- 
jectory. For all the calculations stellar radii R = 10 6 cm 
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Figure 1. Comparison of the three regimes for a neutron star 
with a one-second spin period and a magnetic dipole moment of 
10 31 G cm 3 , pointing in a direction 60° from the line of sight. The 
solid line represents the model containing both QED and plasma- 
based effects, the dot-dashed line shows integrations conducted 
using only plasma effects, and the short-dashed lines represents 
integrations conducted with only vacuum-based effects. 



and mass M — 1.4 Mq were considered. In addition, we 
assumed that the NS's magnetic field was dipola r. The out- 
put of single integration is similar to Figure 3 of IHevl et alJ 
(2003). The si and S2 polarizations change considerably as 
the photon travels away from the NS. At approximately 30 
stellar radii, the polarization decouples. The significance of 
the decoupling radius is discussed below. 



3 RESULTS 

3.1 Net Polarization 

The net polarization of the radiation emitted from the NS 
was determined by integrating over many trajectories on the 
neutron star image. We sampled the image using regions of 
equal solid angle using concentric rings about the centre of 
the NS image. This sampling method could be performed 
due to the inherent symmetry in the circular NS image. By 
symmetry, the S2 and S3 components of the polarization are 
zero. Figure shows the frequency-dependent polarization 
signature for a NS with a dipole moment of /i = 10 31 G cm 3 . 
The full-model polarization (solid line) is plasma dominated 
at low frequencies and vacuum (QED) dominated at higher 
frequencies. In the transitory region, the total polarization 
suffers some destructive interference between the two pro- 
cesses - the total polarization is less than either the vacuum 
or plasma-dominated cases. 

The location of the minimum in the transitory region 
can be understood to be the frequency at which the photon 
decouples from the plasma and vacu um effects at the same 
time. Using the adiabatic condition, IHevl fc Shavivl |2C 
determined the criteria for decoupling, 
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Using this expression, the distance at which decoupling 
occurs can be solved for in both plasma-dominant and QED- 
dominant regimes. To first order these expressions are: 
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Figure 2. Polarization signatures for various magnetic field 
strengths, with period of 1.0s. The plot shows polarization signa- 
tures for increasing magnetic fields, with the lowest line showing 
a field with a dipole moment of 2 X 10 31 G cm 3 , each subsequent 
line representing a factor of two increase in the magnetic field to 
the last line, which shows the polarization signature for a mag- 
netic dipole moment of 6.4 X 10 32 G cm 3 Above each curve is the 
magnetic dipole moment in units of 10 31 G cm 3 . In each case the 
dipole points 30° from the line of sight. 



Here /(a) and g(a) denote the functional dependence 
of the expressions on the inclination of the star. As the res- 
onance occurs when the photon decouples from both the 
vacuum and plasma dominated effects at the same time, the 
critical frequency ui c and the critical radius r c at which this 
resonance occurs are easily calculated: 



oc — r x h(a) 
Pi 



r c oc — r x j (a) 

P7 



(12) 



The results of the numerical model were found to be in agree- 
ment with these relations. 

The critical polarization P c (the polarization at reso- 
nance) is related to the to the radius at which photons de- 
couple: 1 — P c oc (_R/r c ) 7 where < 7 < 2. 

Figure |21 shows how the polarization signature changes 
as the magnetic field strength is varied. It is observed that 
as the magnetic field strength is increased, the strength of 
the resonance feature decreases, but the critical frequency 
doesn't change significantly. Varying the rotational period 
of the NS does change the polarization signature. The vari- 
ation is due to the effects of increased plasma density in the 
magnetosphere. The results seen in Figure [3] show that as 
the rotational period is changed the frequency of the dom- 
inant features are changed, but the overall structure of the 
signature remains similar. 

It is possible that th e charge density of the p lasma 
has been underestimated. lOoldreich fc Julian! JTqIjS) esti- 
mate only the net charge density. However, it is possible 
that the total charge density is larger, as electron-positron 
pairs m ay exist in the m agnetosphere. Extending the equa- 
tions of iMezarosl lll992l) into the formalism presented here, 
the total charge density can be increased arbitrarily. Figure 
Elshows that the introduction of additional charge shifts the 
resonance to higher frequencies, and decreases the depth of 
the feature. Again, the overall shape of the feature remains 
similar. Note that increasing the charge density is compa- 
rable to increasing the rotational rate of the NS. At these 
frequencies increasing the total charge density is nearly co- 
variant with increasing the spin rate of the star, but the 
ellipticity of the modes depends only on the net charge den- 
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Figure 3. Polarization signatures for a dipole moment of 10 31 G 
cm 3 and various periods. The upper panel shows the polarization 
signatures for various NS periods. The lines show, starting from 
the most significant minimum on the left periods of 70s, 32s, 
20s,10s, 7s, 3.2s, 2s, Is, 0.32s, 0.2s, 0.1s, and 0.07s. Some of the 
minima are labeled with the period. The lower panel gives the 
results for various total charge densities. The lowest line shows the 
case in which the total charge density is equal to the net charge 
density. Each successive line above represents an introduction of 
an additional component of the total charge density proportional 
to the net charge density. The second lowest line represents an 
25% excess of charge, and each line above that increases the excess 
by a factor of two. The minima are labeled with the ratio of the 
total charge to the net charge. In each case the dipole points 60° 
from the line of sight. 



sity so the results for higher net charge density (increasing 
O.) and higher total charge density are slightly different. 



3.2 Phase Integrations 

It is currently impractical to measure the instantaneous po- 
larization of an object as faint as a neutron star; therefore it 
is important to calculate the polarization signature from NS 
averaged over many spin periods. To do this the polari zation 
must be calculated in a single basis. iHevl et "ail (l2003t) show 
how the polarization can be calculated. Define the latitude 
of the optical axis (observer) above rotational equator is as 
in and the angle between the magnetic axis and rotational 
axis as 7. If the observer's polarimeter is aligned so that 
so,i is along the rotational equator of the NS and so, 2 is 
denned 45° to it, the components of the polarization can be 
determined from the calculations in the other bases 14 



so,i 



Si 



2 (sin 7 sin < 



1 — (cos 7 sin in + cosir sin 7 cos <f>) 
2 sin 7 sin <j> (cos 7 cos in — cos <j) sin 7 sin in) 



(13) 

so, 2 = si — — — : : : — -j — - (14) 

1 — (cos 7 sin in + cos in sin 7 cos <p) 

The angle a as defined in the previous section can be cal- 
culated using the law of cosines. Figure [I] shows how the 
polarization from a star evolves as the phase varies. The 
quadrature sum of the two terms adds to make s to t, as ex- 
pected. 

The phase-averaged polarization can be found by com- 
puting the mean values in so,i and so, 2 with the average 



14 The expression for so, 2 in lHevl et alj [2003) is incorrect by a 
factor of sinijj. Note that the quadrature sum of so,i an d s O 2 
should be unity. This result is corrected here. 
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Figure 4. Polarization from Neutron Star as a function of phase 
for a NS with 7 = 30° and i R = 30°, for w = 10 13 rad/s and 
fjj = 10 31 G cm 3 . s tot is the total polarization at the particular 
phase, sq 1 is the component of polarization measured in direc- 
tion perpendicular and parallel to the rotational equator of the 
NS, and so, 2 is the component 45° to the equator. 
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Figure 5. Phase-integrated polarization signatures. The upper 
panel shows polarization signatures for various in with 7 = 30° . 
The lines show how the signature changes from sin iji = (bottom 
line) to sinifl = 1.0 (top line, indistinguishable from the x-axis) 
with each above representing an increase in sini^ of 0.1. The 
curves are labelled with the value of sini^j. The lower panel shows 
the results for various 7 with ijj = 30°. The lines show how the 
signature changes from cos 7 = 1.0 (bottom line) to cos 7 = (top 
line) with each above representing an decrement in cos 7 of 0.1. 
The curves are labelled with the value of cos 7. In both panels 
a dipole moment of 10 31 G cm 3 and an NS period of 1.0s were 
assumed. 



value of so, 2 being zero by symmetry. The lower panel of 
Figure |5] shows how the signature changes as in varies. Ob- 
serve that when sin in = 1 (top line) - the observer is look- 
ing directly down on the rotation axis - there is no observed 
polarization. However, for observations at lower latitudes 
there is a significant net polarization and a resonant feature 
at approximately 3 x 10 13 Hz. Statistically, half of all NS 
will have inclination angles such that sini_R ^ 0.5. All the 
models presented in this range show strong features in this 
range considered. 

The upper panel of Figure shows how the signature 
changes as 7 varies. The strongest resonant feature occurs 
when the magnetic and rotational axes are aligned (7 = 
0). However, strong resonant features can be observed for 
cos 7 ^ 0.5. 



4 DISCUSSION 

Of particular attention is the wavelength at which the res- 
onant feature occurs - the near infrared. As this energy, it 
would be possible to measure the polarization of the thermal 
radiation from a NS using a large ground-based telescope 
or Hubble's NICMOS instrument. This instrument is only 
sensitive to wavelengths less than 2.1 /im, so it is possible 
that a large portion of the behaviour could not be probed. 
If the total plasma density is significantly larger than the 
Goldreich- Julian plasma density, the feature moves into near 
infrared or optical energies, and more observing options be- 
come available. 

This letter presents a model for the evolution of po- 
larized light through the magnetosphere of highly magne- 
tized neutron stars. We found that the polarization has a 
frequency-dependent signature, where the total polarization 
decreases significantly in the optical or near IR, dependent 
on the plasma density near the neutron star and the struc- 
ture of the magnetic field. Observations of neutron stars in 
this frequency region may provide useful probes of the prop- 
erties of neutron-star magnetospheres and elucidate the pro- 
cesses responsible for pulsar emission. 
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